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We present a quantization procedure for the electromagnetic field in a circular cylindrical 
cavity with perfectly conducting walls, which is based on the decomposition of the field. A 
new decomposition procedure is proposed; all vector mode functions satisfying the boundary 
conditions are obtained with the help of this decomposition. After expanding the quantized 
field in terms of the vector mode functions, it is possible to derive the Hamiltonian for this 
quantized system. 
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1 Introduction 



The behavior of an atom depends on the mode structure of the electromagnetic field sur- 
rounding it The mode structure in a cavity is determined by its boundary conditions 
which are quite different from the periodic boundary conditions used in Lorentz-covariant 
electrodynamics. There are many interesting effects on the interaction of the atom with the 
field with various boundary conditions. It is possible to calculate them and compare them 
with experimental data 0. Thanks to many new "high-tech" experimental instruments, 
cavity electrodynamics has become one of the most active branches of physics in recent 
years. 

In quantum electrodynamics, we use the plane waves satisfying the periodic boundary 
condition whose quantization volume eventually becomes infinite. As a result, the orientation 
of the field is not essential. On the other hand, the field in a cavity must satisfy the boundary 
conditions and the orientation of the field with respect to cavity walls becomes an important 
factor. For this reason, we have to use vector mode functions to deal with this problem. 
These functions satisfy the boundary conditions and the boundary effect in turn is contained 
in the functions. Thus, for a given cavity, the most important step in the quantization process 
is to construct the vector mode functions. 

Indeed, several authors have investigated the quantization procedure for the field in 
cavities 0]. They have studied the field quantization, the emission rate, and the atomic level 
shifts in the cavities with one or two infinite plane mirrors. Cavity electrodynamics with 
rectangular coordinate system does not present further significant mathematical problems. 
However, the case is quite different for cavities with curved boundary conditions. Even 
for circular cylindrical or spherical cavities, the quantization procedure has not yet been 
firmly established. The purpose of the present paper is to discuss in detail the mathematical 
problems in constructing the vector mode functions with a circular cylindrical boundary 
condition. We hope to discuss the spherical cavity in a later paper. 

Recently, we have carried out the field quantization for several different rectangular 
boundary conditions including a rectangular tube using the vector mode functions 0]. These 
functions have been derived with the help of an orthogonal matrix. However, the procedure 
developed there is not applicable to other cavities in a straightforward manner. 

In this paper, we present a more general method of obtaining the vector mode functions 
and apply it to the field quantization in a circular cylindrical cavity with perfectly conducting 
walls. This will be summarized in two theorems. 

With this in mind, we first decompose the field into three parts in the circular cylindrical 
coordinates in the following section. Then, in Sec. ^, we apply it to the circular cylindrical 
cavity. After obtaining all vector mode functions by using the decomposition, we arrive at 
the quantized field and the quantized electromagnetic Hamiltonian in Sec. ^ 

2 Decomposition of Electromagnetic fields 

Let us now derive the decomposition formula for the electromagnetic field from Maxwell's 
equations. Using the cylindrical coordinate system, we shall decompose the solutions of 
Maxwell's equations into three terms, in preparation for the field quantization within the 
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circular cylindrical cavity in Sec. ^ We emphasize here again that this decomposition process 
is the key to the quantization. 

Maxwell's equations for the electric field E and the magnetic field B in free space are 
given by 

V-E = 0, (2.1) 
V-B = 0, (2.2) 

V xB + dtB = 0, (2.3) 

V X B - \dtB = 0, (2.4) 

where c is the velocity of light in free space and dt = d/dt. As is well known, it follows from 
Maxwell's equations that the fields E and B satisfy the wave equations: 

A-19,^)e = 0, (a-19,^)b = 0, (2.5) 

where A is the Laplacian operator. 

The electromagnetic field E and B can be written in the cylindrical coordinates (r, z) 

as 

E = Et + E„ B = Bt + B^, (2.6) 

where E^ = GrEr + e^E'^ is the transverse component of the field and E^ = e^i?^. Here 
e^, e<^, and are the unit vectors in the r, ip, and z directions, respectively. Similarly, 

V X E, V ■ E, and V0 (0 a scaler function) can also be divided into two parts. For example, 

V X E = Vr X E + X E, where the first term contains the derivatives with respect to r 
and V9, while the second term contains the derivative with respect to z. For simplicity, the 
derivatives with respect to r, y?, and z are described as d/dr = dr, d/d^ = and d/dz = dz, 
respectively. 

Since x E^ = 0, Eq. (O) gives 



dtBz = -Vt X Et, dtBr = -Vt x E, - V, x Et. (2.7) 
Similarly, from Eq. (|2.4|), we have 



^ dtEz = Vt X Bt, \dtET = Vt x B, + V, x Bt. (2.^ 



Equations (Ij) and (|]|) give 



c 



dfET = Vt X dtBz - V, x (Vt x E, + V, x Et) 



^O^Bt = -^VtxS^E, -V, X (VtxB, + V. xBt), (2.9) 



which leads to 



— 9?Et + V^ X V^ X Et - Vt X Vt X E^ = V X dtB^ - V x V x E^, 

^^^Bt + V. X V, X Bt - Vt X Vt X B, = --^V x S^E, - V x V x B„ (2.10) 



where we have used Vt x = V x E^. Let us note that 

V, X V, X Et = -A,Et, 
Vt X Vt X E^ = -AtE^, 
1 



— 9;Et = (AE)t = AEt, (2.11^ 



where (AE)t is the transverse part of AE; the last equation holds in the present coordi- 
nates, although it is not correct in general. (For example, it is not correct in the spherical 
coordinates.) Then, from Eq. (p. 101) , we arrive at 

AtE = -VxVxE^ + Vx dtB,, 

AtB = -4v X dt^, - V X V X B^. (2.12) 

To rewrite Eq. ( |2.12|) , we must obtain an expression for the components Ez and B^. 
Suppose that the field is in a finite region and expand and in terms of a certain 
complete system of functions with mode s: 

E,{r,t) = ^(E,,(r,t) +C.C.), 

s 

B,{r,t) = ^(S,,(r,t)+c.c.), (2.13) 

s 

where 

E^r, t) = E,,(r)e-^"-*, 5,,(r, t) = B,,{v)e-'-^^\ (2.14) 

Here Uga ijjJsa > 0, a = 1, 2) is determined by using given boundary conditions. Since E-^ 
and Bz satisfy the wave equation, their components satisfy the Helmholtz equations: 

/\Ezs = -kl^Ezs, ABzs = -kl^Bzs, (2.15) 

where 

kl = (2-16) 

We assume that the components satisfy 

dlEz, = -hl.Ezs, dlBz, = -hl^Bzs, (2.17) 

where h^^ is also determined by the boundary conditions. Then we have two dimensional 
Helmholtz equations: 

^tEzs = -gsiEzs, ^tBzs = -gs2Bzs, (2-18) 

where 

gl = kl-hl. (2.19) 
Now we define two functions E„ from E^s and B^s with g"^^ 7^ as 

F^{r,t) = J2 [Fsa{r,t) + c.c.] = ^ [F,.(r) e-^^^^* + c.c], (2.20) 
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where 



Fsi = Ezs/gl^, Fs2 = Bzslgl2- 



[2.21] 



The functions F^. and their components Fga satisfy the same equations as the z components 
of the field: 



1^ 



AtF, 



-0^ F 



(2.22) 



From Eqs. ( 2.21 ) and ( 2.22| ) the component Fga is a solution of the Poisson equation 



and then the functions At-Fo- satisfy 



(2.23) 



ArFi = -Y. {g^Fsi + c.c.) = - ^ (F,, + c.c.) 
AtF2 = - Y1 i9^s2Fs2 + c.c.) = - J2 (Bzs + C.C.; 



(2.24) 



On the other hand, if there is a component E^s or B^g with g"^^ = 0, Eq. ( |2.18| ) reduces 
to the two dimensional Laplace equation: AtE^s = or AtB^s = 0. Now define Eq^ and 
Bqz as 

Eoz = + c-c), = (Bzs + c.c), (2.25) 

which satisfy 

AtEoz = 0, ArSo, = 0. (2.26) 

Then we have a useful expression for E^ and B^; for any F^ and B^ there exist functions 
Fa- given by Eq. (|2.20| ) and functions Eq^ and i^oz are given by Eq. (|2.25|) such that 



F. 



-ArFi + F, 



02, 



5. 



/\tF2 + -Bi 



02- 



(2.27) 



It is worth emphasizing again that the functions F^ are constructed with the components 
Ezs and Bzs with g^^^ ^ 0, while Eqz and -Bq^ consist of the components with gl^ = 0. 
Using Eq. ( p.27|) and defining Fo- by 



e F 



we can rewrite Eq. (p.l2|) as 



At(E - VxVxFi + Vx dtF2) = -V x V x Eq^ + V x dtB. 



At ( B 



1 



V X dtFi - V X V X Fa 



1 



V X aEn^ - V X V X B 



0X5 



(2.28) 

(2.29) 
(2.30) 



where Eqz = GzEoz and Bqz = e.zBoz- Let us define Eq and Bq as the quantities in the 
parentheses at the left hand side in Eqs. ( p.29|) and ( p.30|) , respectively. The results of this 
section are then summarized in the following theorem. 
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Theorem 1: The field can be decomposed into three components as fohows: 

E = VxVxFi-Vx dtF2 + Eq, 



B = X dtFi + V X V X F2 + Bo, (2.31) 



c2 



where Eq and Bq satisfy 



AtEo = -V X V X Eo. + V X dtBo,, 



ArBo = -^V X dtEo, - V x V x Bq.. (2.32) 



C2 



Theorem 1 plays a central role in performing field quantization in this paper. Take the z 
component of Eq. (p.31| ), we find that the z components of Eq and Bq are given by Eq^ and 
Bqz in Eq. (|2.27|) , respectively. Furthermore, Eq. ( p.32|) is consistent with Eq. ( p.l2|) . 

Let us consider the physical meaning of the above decomposition formula. First it is 
worth emphasizing that each term in Eq. ( p.31|) is a solution to Maxwell's equations. That 
is, each term containing F^r {F^ term) satisfies Eqs. ( p^.l| ) - (O), which is easily shown by 



using the fact that F^ satisfies the wave equation ( |2.22|) . As a result, the third term Eq, Bq 
is also a solution to Maxwell's equations, so that it satisfies Eq. ( p.l2[ ), as mentioned above. 

Next take a particular example of the boundary condition satisfying Eqz-, Bqz = 0. Then 
the field Eo,Bo describes so-called the TEM (transverse electromagnetic) mode p. Since 
the Fi term has no z component of the magnetic field, it becomes TM (transverse mag- 
netic), while the F2 term becomes TE (transverse electric), because it does not contain the 
z component of the electric field. 

In the above particular case where Eq^, -Bqz = 0, Maxwell's equations for the field Eq, Bq 
reduce to 

Vt ■ Eo = 0, Vt X Eo = 0, 
Vr ■ Bo = 0, Vt X Bo = 0, 
V, X Eo + dtBo = 0, 

V. X Bo - \dtBo = 0. (2.33) 



3 Determination of Functions F^j 

We are considering here the cavity enclosed by a circular cylindrical wall with radius a and 
height L: r < a, < z < L. We assume that the cavity has perfectly conducting walls 
at 2; = 0, L and at r = a. The tangential component of the electric field E|tan and the 
normal component of the magnetic field B|norm must accordingly vanish at the boundary of 
the cavity. 

The above boundary condition reduces to that for the z components 

Ez = 0, drBz = 0, (r = a), (3.1) 
B, = 0, d,E, = 0, {z = 0, L). (3.2) 
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It is easy to get the second condition in Eq. ( |3.1D if we take the component of Eq. 



Equation leads to the second condition in Eq. ( p.2|) 

Before obtaining the components E^^ and B^s-, we show that the folowing lemma. 

Lemma 1: Eq = 0, Bq = 0. 

Proof: First we show that Eqz,Bqz = 0. Consider the component Eq^, which satisfies 
the Laplace equation AtEq^ = [see Eq. ( p.26|) ] with the boundary condition Eq^ = at 



the boundary (r = a). Then we have Eq^ = 0; this is a well known property of the Laplace 
equation. 

Next consider the component Boz satisfying AtBqz = with the boundary condition 
drBoz = (r = a). Then Boz must be a constant, i.e., independent of the variables r and ip: 
Bqz = f{z, t). Applying the two-dimensional Gauss' theorem to V-Bq = Vr-Bor+f^^-Boz = 0, 
we get OzBqz = 0, because 5or = (r = a) and, as a result, Bq^ is independent of the variables 
r, ip, and z. On the other hand, Bq^ must be zero at ^ = 0, L, which leads to Bq^ = 0. 

As mentioned in the preceding section, the field Eq, Bq becomes TEM satisfying Eq. 
(|2.33| ), because Eq^, Bq^ = 0. There exists no TEM in the present cavity. That is, from Eq. 



( |2.33| ) we have Vt x Bq = and Vt ■ Bq = 0. As a result, there exists a function ip such 



that Bq = — VtV^, which gives Bq = by the boundary condition at r = a. Similarly, we 
get a vector function A such that Eq = V x A. The boundary condition leads us to Eq = 0. 
Q.E.D. 

Since Eq, Bq = 0, Theorem 1 shows that the field is constructed from the functions 
Fo-. Consequently, from Eq. ( |2.20| ), the field only contains the components with g"^^ ^ 0. 
Moreover, it is easy to prove that g1^ > by using Gauss' theorem and the boundary 
condition. 

Next we solve the Helmholtz equation ( p.lSj) for the components E^s and B^s under the 



boundary conditions ( ^.l] ) and (^.21). Since g^^ > 0, the solution is given by 

EUr,t) = Csi{t)JUXsir/a) d"^^ cos{mi z / L) , 

Bzsir, t) = Cs2it)Jrnixs2r/a) e'^^ sin(n7r^/L), (3.3) 

where the mode index is s = (m, /i, n) {m = 0, ±1, ±2, ■ ■ ■ ; fi = 1, 2, 3, ■ ■ ■ ; n = 0, 1, 2, ■ ■ ■), 
Jm is the Bessel function of the first kind, Xsi = Xm^ti is the /xth zero point of Jm, and 
Xs2 = XmiJ.2 is the fith zero point of J^, the derivative of Jm- It follows from Eq. ( |2.14D that 

From the solution ( ^^) we have 

AtEzs = -{kl^ - {n'K/Lf)Ezs = -{xsi/afEzs, 

AtBzs = -{kl^ - {mx/Lf)Bzs = -{Xs2/afBzs, (3.4) 

which gives 

gl = iXsa/af, kl = {xsjaf + {riT^/Lf. (3.5) 

Let us next obtain the functions F^-, which are defined in Eqs. ( |2.20|) and ( |2.21| ); they 
are given by 



F,i(r,t) = 2 = ^\h^ — asi[t)4)si{j) 

Yd V ^Sn 
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s2 



Xs2 



2^0 



as2(t)^s2(r) 



(3.6) 



where we have introduced a^^ and ipsa'- 

as.{t) = a,,(0)e-^-='^*, (3.7) 
^psi{r,t) = CsiJm{Xsir/a) e'"''^ cos{nnz/L), (3.8) 
^.2(r,t) = c,2JmiXs2r/a)e'"'^sm{niTz/L), (3.9) 

where Csa are normahzation constants. The functions ipsa have the orthonormahty property, 
which is used in quantization in the next section. 

Lemma 2: 



1, 



(3.10) 



where dr = /cavity rdrdipdz, V is the cavity volume, and 

«s2 = Oim^i2 = JmiXm^i2) — Jm+liXm^i2) ■ (3-11) 

Here the quantity cos{nnz/L) in Eq. ( |3.8| ) must be changed to 1/V2 when n = 0. 

Proof: The orthonormahty property (|3.1CI| ) is easily derived from the following equality 
for the Bessel functions: 



1 

drrJm{Xmfiar/a)Jm{Xmf,'ar/a) = -a'^asaSf.f,'. 

To prove Eq. ( p.l2| ), let us first observe the following two integrals 

1 



(3.12) 



V? — V? 

A2 A. J 



At Aj 



XiJm+l{Xi)Jm{Xj) - XjJm+l{Xj)Jm{X'. 
~ XiJrniXi) JmiXj) + Xj JmiXi) JmiXj) 



(3.13) 



dx xJ,^(^Xix) 



1 r 
2 

1 r 



'^miXi) -Jm—liXi) '^m+l{Xi 



JnXXi) Jm+l{Xi) J rnkXi) ^ m+liXi) 



(3.14) 



where Xi (^ = 2, ■ ■ ■) is a positive real number and Xi Xj- 

Consideri ng tha t Xm^ i is th e /xth zero point of Jm{x) and that Jm-iiXm/ii) + Jm+i{Xm,ii] 
0, from Eqs. ( |3.13D and (|3.14| ), we have 



1 1 

dx xJfYii^Xm^l-^^ Jmi^Xmfi'l-^^ ~OirnpLl ^ i 



(3.15) 



where Om^i is given by Eq. ( |3.11|) . If we change the variable x to r/a, Eq. ( p.l5|) gives Eq. 
(|3:T^ ) for a = l. 

Similarly, Eqs. ( p.l3|) and (|3.14|) lead to 



1 1 

dx xJ,^(^Xm,fi2-^) 'JiniXmiJ,'2-^) ~0irn^2 ^ ^fj,' i 



(3.16) 



where ^^^2 is given by Eq. ( p.ll| ). Equation ( p.l6| ) gives us Eq. ( |3.12| ) for a = 2. Q.E.D. 
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4 Vector Mode Functions and Field Quantization 

After defining the vector mode functions, we quantize the field and obtain the quantized 
Hamiltonian. The mode functions are constructed according to the decomposition formula 
given in Eq. ( p.31| ) 

In the cavity with the perfectly conducting walls, the decomposition (|2.31| ) in Theorem 
1 reduces to 



E = VxVxFi-Vx dtF2, 
B = — V X dtFi + V X V X F2, 

because of Eq, Bq = [see Lemma 1]. 

Substituting the functions F^. in Eq. ( ^.61 ) into E in Eq. ([4.1D, we find 



(4.1) 



E(r,t)=^E 



2eo 



a,^(t)u,^(r) - a*^(t)u*^(r) 



where the vector mode functions u^o- are given by 

Usl = V X V X e^lpsl, Us2 = iuJs2^ X ez'^s2- 

Similarly, substituting into B and considering the equahty 

V X u,i = kl^V X e^?/',i, 

we obtain 



B(r,t) = E 



asa{t)V X Us^(r) + a*^(t)V x u*^(r) 



(4.2) 

(4.3) 
(4.4) 
(4.5) 



Before obtaining the quantized field and Hamiltonian, let us present some properties of 
the mode functions. 



Lemma 3: 



0, 



:a + e)u,, = 0, 



u 



sa I tan " 7 



0, V X u 



S(T norm 



0, (on walls). 



rfr u*^(r) ■ uy^/(r) = 5ss' S^a', 



(4.6) 
(4.7) 
(4.8) 

(4.9) 



Proof: We give a proof of the orthonormality property ( |4.9| ), because it is obvious to 
prove the other equations. From Eq. (|4.3|), the mode function u^i is rewritten as 



(4.10) 



which gives 



^ll^ll - KlKl flits' 



2 u2 



+ V 



s'l 



(4.11) 
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Taking into account Gauss' theorem and the boundary conditions, we find that the last two 
terms in Eq. ( |4.11| ) have no effect on the volume integration of u*^ ■ Us'i. Equation ([4.11|) 
gives the orthonormality relation for u^i: 



where we have used Eq. (|3.1CI| ) and set 



(4.12) 



In the case of 0^2, we make use of the equality 

(V X e^i^l^) ■ (V X e,ips'2) = e.tpl^ ■ (V x V x e^^./s) + V ■ le.tpl^ x (V x e,^,/2 



(4.13) 



Xs'2 



(4.14) 



The orthonormality property of is then obtained from Eq. ( ^.lUD if we set the normal- 
ization constant: 

I ^ 

V ^as2Xs2^s2 

Finally we show that the mode functions u^i and Us'2 are orthogonal to each other. Since 
has no z component, we see that 

u*i ■ u,/2 = iujs'2[kl^e^ipl^ + Vd^ipl^^ ■ (V x e^iljs'2) 

= iujs'2'^{dzip*si) ■ (V X e^ilJs'2) 

= iuJs'2''^ ■ [{dzip*si){V X ezilJs'2)], (4.16) 

which has no effect on the volume integration of u*-^ ■ Us'2- Q.E.D. 

We are now ready to carry out the field quantization. To get the quantized field the 
functions asa{t) are regarded as annihilation operators, which annihilate photons with indices 
sa and satisfy the commutation relation: 

[asa{t),al,^,{t)] = 6ss'6aa', (4.17) 

where al,^,(t) are creation operators. Then the field given by Eqs. ( |4.2| ) and ([4.5|) becomes 
operators. The quantized electromagnetic Hamiltonian is obtained by using the equality 

^dr(V X u:,(r)) ■ (V x u,,.,(r)) = ^rfr u:,(r) • u,,.,(r). (4.18) 

These are summarized in the following theorem. 

Theorem 2: The quantized field and the Hamiltonian are given by 



E(r,t)=2E 

sa 



2en 



asa{t)^sa{r) - aj„(t)u*^(r) 



(t) V X u^„(r) + aL(t)V x u*^(r) 



4. 



(4.19) 

(4.20) 
(4.21) 
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Theorem 2 will play a fundamental role when we investigate the electromagnetic property 
of a physical system in the circular cylindrical cavity. There are also important issues in 
the present cavity including, for example, spontaneous (stimulated) emission, absorption, 
the atomic energy shifts, and the Casimir effect; they will be investigated with the help of 
Theorem 2. 

It is necessary to expand the field in terms of the vector mode functions as in Theorem 
2. According to Fourier analysis, the field can be expanded in terms of, for instance, the 
periodic functions as in ordinary quantum electrodynamics. However, in this case, the 
operators a^o- and al^ cannot describe the physical photons. Consider the situation where 
the field consists of only one mode corresponding to the photon with aga'^sa, where Uga 
satisfies the periodic boundary condition. This must be impossible, because this field does 
not satisfy the boundary conditions; the photon is then fictitious. 



5 Conclusions 

The quantization procedure of the field in the cavity without Eq, Bq is as follows: (a) ob- 
tain the decomposition formula ( p.31| ) in the circular cylindrical coordinates; (b) solve the 



Helmholtz equations ( p.l5|) for the components Ezs and Bzs under the boundary conditions 



( |3.1| ) and ( |3.2| ); (c) determine the functions Fo-, which are the solutions of the Poisson equa- 



tion (|2.23|) , from E^s, Bzs] (d) substitute F^- into the decomposition formula (|2.31|) and obtain 



the vector mode functions satisfying the orthonormality property ( [4.9| ); (e) then we arrive 
at the quantized field and Hamiltonian in Theorem 2. 

As has been shown explicitly in Theorem 2 in Sec. 0, the energy hug^^ of the photon 
depends on the polarization index as well as the wave number index. This is one of the 
characteristics of the circular cylindrical cavity. The photon energy is independent of the 
polarization index in a rectangular cavity and in free space. 

In the entire process of quantization, the decomposition formula in Theorem 1 plays the 
essential role. We must find such a formula in an appropriate coordinate system if we study 



other types of cavities. Also, if a cavity allows the third term Eo,Bo in (|2.31|) , we have to 
construct a methodology to deal with the problem. We hope to discuss the quantization of 
electromagnetic fields inside a spherical cavity in a future publication. 
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